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Abstract. In recent years, theories in which the Einstien-Hilbert lagranigan is 
replaced by a function f{R) of the Ricci Scalar have been extensively studied in 
four-dimensional spacetime. In this work we carry out an analysis of such theories 
in two-dimensional spacetime with focus on cosmological implications. Solutions 
to the cosmological field equations are obtained and their properties are analysed. 
Inflationary solutions are also obtained and discussed. Quantization is then car- 
ried out, the Wheeler-DeWitt equation is set up and its exact solutions obtained. 
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1 Introduction 



Attempts to modify the theory of general relativity, by including higher order 
invariants in the action, started not too long after its inception [1], [2]. Later the 
non-renormalizability of general relativity gave impetus to the inclusion of higher 
order terms in the action [3]. More recently it was shown that when quantum 
corrections are taken into consideration, higher order curvature invariants need to 
be added to the low energy gravitational action [4]. Such considerations further 
increased the interest in constructing theories in which the Einstien-Hilbert action 
is extended by the inclusion of higher order curvature invariants with respect to 
the Ricci Scalar. Our interest here is in the so-called f{R) theories of gravity. In 
these theories the Lagrangian in the Einstein-Hilbert action 



where n = SttG, G is the gravitational constant, g is the determinant of the metric 
tensor and R is the Ricci scalar (in units c = h= 1), is generalized to become 



In eq.© f{R) is a general function of R [5]. Our focus here is on the cosmologi- 
cal aspects of f{R) theories. 

Now in another direction, the quest for quantum theory of gravity has led to the 
study of the simpler case of gravitational theory in two-dimensional spacetime. 
Such a spacetime provides an interesting arena in which to explore some funda- 
mental aspects of both classical and quantum gravity. The reduction in the degrees 
of freedom greatly simplifies the analysis of the field equations. This leads to ap- 
preciable understanding of several problems in gravity theory. In two-dimensional 
spacetime, the two-dimensional gravitational constant G2 is dimensionless and 
formally the theory with the bare action is 



where = 8716*2. is power counting renormalizable in perturbation theory. How- 
ever the Einstein— Hilbert action term is purely topological in two dimensions. In 
fact in two spactime dimensions, the curvature tensor R^uXp has only one indepen- 
dent component since all nonzero components may be obtained by symmetry from 
i?oioi • Equivalently the curvature tensor may be written in terms of the curvature 
scalar [6] 




(1) 




(2) 




(3) 



RfiuXp — 7;R i9fj.x9up — gfipQux) , 



(4) 
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so that R alone completely characterizes the local geometry. Eq.© implies that 



(5) 



so that the Einstein tensor = R^i, - ^g^j,R, vanishes identically and the usual 
Einstein equations are meaningless in two dimensions. This led to various models 
for gravity in two-dimensional spacetime being proposed [7]. Of special interest 
are those models that involve a scalar field, the dilation, in the action [7] — [10]. 
We have previously studied some aspects of classical and quantum cosmology 
in two-dimensional dilation gravity models [11], [12]. In the present work we 
study theories as an alternative way to formulate gravitational theory in two- 
dimensional spactime and explore some of their cosmological implications. We 
restrict ourselves to the classical theory only. 

In sect. 2 we set up the f(R) gravity theory in two-dimensional spacetime and 
derive the general field equations. We then specialize to the case of the Friedmann 
— Robertson —Walker metric and obtain the field equations with matter treated as 
a perfect fluid. Sect. 3 is devoted to obtaining solutions to the cosmological field 
equations under various conditions of matter or radiation dominance. Properties 
of these solutions are discussed in sect. 4. In particular, conditions for ensuring 
cosmic acceleration and solving the horizon problem are elucidated. Inflation is 
discussed in sect. 5 and solutions to the field equations in the absence of matter or 
radiation are obtained and their properties are discussed. In sect.6 we carry out the 
quantization. We establish the Wheeler-DeWitt equation and obtain its solutions. 
In sect. 7 we offer some concluding remarks. 

2 Field equations 

We write the two-dimensional action for f{R) gravity as 



is the gravitational action and Im is the matter action [13]. The field equations 
can be derived by varying the action with respect to the metric tensor g^,y. Upon 
noting that the stress-energy tensor is defined by 



I = Ig + Im, 



(6) 



where 




(7) 




(8) 



we derive the following field equation 
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f{R) R,, - \g,uf{R) - g^uOfiR) + V^V,/'(i?) = -g,T,, . (9) 

In eq.dH) R^^ is the Ricci tensor, the prime denotes the differentiation with respect 
to R and the operator □ is defined by 

nf{R) = -1= d,iV^ g'^'' dJ'iR)) . (10) 
Using eq.© we can write eq.® as 

^g,,{f'{R)R - f{R)) - g,.af'{R) + V,VJ\R) = -g,T,, . (11) 

In the following we shall be concerned with cosmological implicatons of eq.lfTTI). 
For this purpose we shall adopt the Friedman-Robertson-Walker (FRW) metric 
which in two-dimensional spacetime reads ( c = 1) 

= -dt^ + dx' , (12) 

1 — k 

in terms of the comoving coordinates x and t. The quantity a{t) is the usual time- 
dependent cosmic scale factor. A change of variable dx'^/ (1 — kx'^) — )■ dx"^ leads 

ds^ = -dt^ + a^{t) dx^ . (13) 

Thus in two dimensions the time evolution of a(t) is not affected by the value of 
/c = 0, ±1 corresponding to the three different cosmological models [14]. This 
is unlike the four-dimensional case. The values A; = 0, — 1 still describe spatially 
open flat and hyperbolic universes respectively while k = 1 describes a closed 
universe. The stress-energy tensor of the homogeneous isotropic universe is taken 
to be that of a perfect fluid 

Tf,u =pgf.u+ {p + p) , (14) 

where p is the pressure, p is the energy density and is the comoving velocity. 
Using eq.(fT3l) and eq. (fT4l) we obtain from eq. (fTTI) the following two independent 
cosmological field equations 

liR f{R) - f{R)) + -dJ'{R) = g,p , (15) 
^{R f\R) - f{R)) + dlf\R) = -g^p, (16) 
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where we use the dot as well as dt to indicate differentiation with respect to time. 
We note that if f{R) is expressed as a sum of powers i?" of R, then a term linear 
in R would cancel out in the bracketed terms in eq.(fT5l) and eq.(fT6l) and would not 
contribute to the derivative terms either. Hence it has no effect on the dynamics. 
The stress-energy tensor obeys the conservation law 

V"T„^ = , (17) 
and this, for a perfect fluid, gives rise to the following two equations 

^7" V„p+ (p + p) V"f/„ = 0, (18) 

and 

(p + p)f/" v«f/;3+ {gap+ u^u^)V''p = o. (19) 

For the FRW metric of eq.(fT3]) one readily obtains from eq.(fT8]) that 

^{pa) = -p. (20) 
da 

Assuming an equation of state of the form p = 7p where 7 is a constant, eq. (|20l) 
immediately leads to 

p = Ca'"'-^ , (21) 

where C is a constant. Eq.(fT9l) is seen to be identically satisfied and does not 
give rise to anything new. For a pressureless (dust) pure matter universe (p^ 7^ 0, 
Pr = 0, 7 = 0) we have 

Pm = Cm a^-^ , (22) 
while for a pure radiation universe (p^ = 0, p,. 7^ 0, 7 = 1) one has 

Pr = Cr , (23) 

Denoting the present time by to and using the usual notation of ao = a(to) and 
po = p{to) to denote present-day values of these quantities, we can write for a 
matter-dominated universe 

Pm = 0, pm{t) = pmO -7TT , (24) 

a{t) 

while for a radiation-dominated universe one has 



Pr = Pr{t) = pro[ ^ ) , (25) 
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Finally we wish to note that for the FRW metric the curvature scalar of this two- 
dimensional universe is given by 

R = -^,, (26) 
a[t) 

where a = ^ . 



3 Solutions of the cosmological field equations 

In this section we seek solutions of the cosmological field eqs.(fT5]) and (fT6l) with 
the energy density and pressure given by eqs. (|24l) and (1251) for each component 
of the cosmological fluid thus obtaining two sets of equations. For the matter 
dominated epoch we obtain the following equations 

i {Rf'iR) - f{R)) + ^dJ'iR) = g^p^o^ , (27) 

^{RnR)-f{R)) + d^nR) = 0. (28) 
For the radiation dominated epoch the corresponding equations read 

I {Rf{R) - f{R)) + -d,f'{R) = g^pro^ , (29) 

^ {RfiR) - f{R)) + dU'{R) = -a.Pro^ , (30) 

To proceed further we need to specify the function f{R). Similar to the proce- 
dure followed in the four-dimensional case [5] we take for f{R) the following 
expression 

f{R) = R + aR"" , (31) 

where the real constants a and n are, at this stage, only restricted by a 7^ and 
n ^ 1. Upon substitution of eq.dSTI) into eqs.dTTl) and (l28l) we obtain the equations 

hn- l)aR^ + n(n - l)aR''~^R- = g^Pmo- , (32) 
^ d d 

2nRR + 2n{n - 2)R? + R^ = . (33) 

Equations (l32l) and (1331) describe the matter dominated epoch and we shall attempt 
to find solutions for them now. We start with eq.(l33]) and note that in terms of the 
function z{R) defined by 
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z{R) = , 

the equation is transformed into the following form 

dz 2(n - 2) 

This equation is easily solved and we obtain for n ^ ^ 



(34) 



(35) 



z{R) = R^ = - 



-R^ + CiR 



4-2n 



n{2n - 1^ 

where Ci is a constant. Eq.(l36l) then leads to the parametric solution 



(36) 



t = ± 



R' 



n{2n - 1) 



CiR 



4-2n 



dR + C2, 



(37) 



where C2 is a constant. For n = 2 and Ci 7^ one can carry out the integration 
using the result [14] 



dx 



X 



(ir-X"+2)V2 



1 a + 3 a;' 



a+2 



2'a + 2'a + 2' isT 



(38) 



where a, K are constants and 2-F1 is the hypergeometric function. We obtain 



(39) 



Ideally one should solve eq.(l39l) to obtain i? as a function of the cosmic time t and 
plug that into eq.(l32l) in order to solve for a{t) in the case of n = 2, but that is 
a difficult task. Instead we consider solutions for which Ci = in eq.(l36l) and a 
general ^ 7^ |. One can then easily derive that 



R = - 



An{2n 



(40) 



where we have renamed the integration constant C2 as tm- In fact one can verify 
directly by substitution that the expression for R in eq. (|40l) is a solution of eq. (|33l) 

Next we substitute eq.(l40l) into eq.(l32l) and obtain 



+ Ai{t- tm) = Kao {t - tm) 



2n-l 



where 



Ai = ^{n - 1)A^"; A2 = -2n{n - 1)A^"-^ 



(41) 
(42) 
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and 

N = An{l-2n)- k = . (43) 

a 



We readily solve eq.(l4Tl) and get 

a{t) = C(t - + K{t - , (44) 

where C* is a constant and 

K = (45) 

Clearly n must be such that Ai and A2 are real and K is finite. We will return to 
this issue later. It is interesting to note that the t dependence of i? is (t — tm)^^ 
and thus independent of n, while that of a(t) does depend on n. We also note that 
the relation R = — ^ is satisfied by the solutions for R and are given in eqs.(l40l) 
and (l44l) respectively. We further note that the second term in eq.(l44l) is a solution 
of eq. (l4T]) in its own right. On the other hand the first term in eq. (l44l) is a solution 
of the homogeneous form of eq. (|4T]) . Furthermore the constants C and K must be 
such that a{t) is positive. 

We now turn to the case of radiation. Upon adding eqs.(|29l) and (l30l) we obtain 
the equation 

Rf'iR) - f{R) + dlfiR) + -dtfiR) = . (46) 

a 

Employing in eq.(l46l) the expression for f{R) given in eq.(l3TI) above yields 

nRR + n{n - 2)R^ + R^ + nRR- = . (47) 



Next we use eq. (pTI) in eq.(l29l) and obtain 

^{n - l)ai?" + n{n - l)ai?"-2i?^ = g^Pro^ . (48) 

Motivated by the structure of the solutions for the cosmological equations in the 
case of pure matter above, we seek solutions for R(t) and a{t) of eqs.(l47l) and 
(|48l) in the form of powers int — tr where U is some reference time. We obtain 
the following results 

R{t) = 2n{n-l){t-try^ , (49) 



a{t) = B{t - trf , (50) 
where the constant B is given by 
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B 



nl/2 



ao 



(51) 



_n{n - 1)(1 - 3n)[2n(l - n)Y-^a 

Note that, as in the case of matter, the t dependence of R{t) is independent of 
n, the only such dependence appears in the overall coefficient. We also note that 
the relation R = — ^ is satisfied by the solutions for R and are given above in 
eqs.(l48l) and (|50l) . For an expanding universe one must have n > \ and B > 
Furthermore the value of n must ensure that the bracketed term in eq. (|5TI) is finite 
and real. 



4 Properties of the solutions 

We now discuss some properties of the solutions of the cosmological field equa- 
tions found in the previous section. Let us first look at the radiation dominated 
case and determine whether our vision of the universe is limited by a particle 
horizon. At a given cosmic time tg the proper distance d(ts) of the emitter is given 
by 

d{t,) = a{Q / — - , (52) 



(53) 



It. a{t') 

Where tg is the time of emission of the photon. Using eg. ([50b we obtain 



d{ts 



n 



We can view tr as signifying the onset of the radiation epoch. We see that as 
te — )■ tr, d(ts) is finite for 1 — n > and diverges for 1 — n < 0. Hence no particle 
horizon problem will arise if n > 1 which is the same condition required for an 
expanding universe. Reality of B also requires n to be an integer. For n an even 
integer, the parameter a must be positive while for n odd, a should be negative. 
Thus we take n to be a positive integer greater than one. Next we note that the 
cosmic acceleration a{t) which is given by 

a{t) = n{n - 1) B {t - t^)""^ , (54) 

is positive for t > tr since n > 1 and is constant for n = 2. Now in two- 
dimensional spacetime the radiation energy density pr oc were T is the tem- 
perature [14] and it follows therefore from eq. (l25l) that 

aocT-^ (55) 

Since we have a — )■ as t — )■ t^, we conclude that this radiation universe has a hot 
big bang origin. 
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Next we turn to the case of the matter dominated universe described by eqs. (I40l) 
and (l44l) . First let us consider the case C = when the scale factor becomes 



a{t) = K{t - t„,) 



2n 



(56) 



As we have stated earlier this is viable because it represents a solution of eq. lHTI) . 
The time tm can be taken to signify the onset of matter dominance. The proper 
distance d(ts) is now given by 



\2n 



1 - 2n 



a-2n 



a-2nl 



(57) 



Hence no particle horizon will arise if 2n > 1. Also as we stated following 
eq.(l45l). the parameter n must be such that the constants Ai and A2 given by 
eq.(l42l) are real. Since for 2n > 1 the number of eq.(l43]) is negative, it follows 
that n has to be a positive integer. Now the requirement that a(t) > for t > 
implies that K > 0. For n even we have ^2 > and hence a should be positive 
to ensure A' > while for n odd one has A2 < and a should be negative. Since 
we exclude n = 1, the smallest permissible value is n = 2. For such values of n 
it is evident that the cosmic acceleration d{t) is positive. Finally we observe that 
for the pure matter universe we have a{t) — )■ as t — )■ t^- 

We now consider the case C ^ 0. Using eq.(l44l) the proper distance is now given 
by 



d(ts) = a{ts) 



2n-l 



\4n-l 



-dt. 



(58) 



C + K{t- tr,^, 

It is clear that the integral converges for tg — ^ tm and we do a have particle horizon. 
Performing the integral we determine the proper distance to the horizon to be 



d(ts 



C 



2n(4n-l) 



ln{l + is 
An - I 



An 



^ 2n-l 
k=l 



2n-K{2k - 1^ 
An - I 



/nil— 2is cos 



2k 

An 



2ri-l 



+ 



An 



sm 



k=l 



'2mr{2k - 1^ 
An - 1 



arctg 



- COS^^TC' 

4n— 1 

An—l 



where 







k) 



(59) 
(60) 
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Let us now study further properties of the solution given in eq.(l44l). In the follow- 
ing we consider only values of t such that t > ■ Now it is evident that, except 
for values of n in the interval < n < |, the first term in eq. (l44l) dominates for 
t near tm when n > ^ while the second term dominates for n < 0. Hence to 
ensure positivity of the scale factor we require that both C and K be positive. For 
< n < ^ , C and K can have opposite signs but only in such a manner so as to 
keep a > 0. We shall for simplicity assume that C > and A' > for all values 
of n. Next we observe that outside the interval < n < ^ , the number of 
eq.(l43l) is negative and to ensure the reality of A2 given by eq. (l42l) . the number n 
has to be an integer. We readily deduce that for a > 0, n can be a positive even 
integer or a negative odd integer. On the other hand for a < 0,n can be a positive 
odd integer or a negative even integer. The cosmic acceleration d{t) is given by 

'd{t) =2n{2n-l){t-t,r,)-^a{t) . (61) 

It is seen that a < for < n < ^, a = for n = ^ and a > for n < 
or n > I . 

Next we consider the behavior of a{t) as t — )■ for the case C ^ 0. We see from 
eq. (l44l) that for < n < |, a{t) — )• as t — )■ tm and accordingly the temperature 
T — !■ 00 in this limit. For n = |, we have a(t) C as t ^ tm and T is finite. 
However for n outside the interval < n < | the behavior of a(t) is very different 
as t — )■ tm- We see that a(t) — )■ 00 in this limit and and energy density pm as well 
as the temperature tend to zero. As t increases beyond the value tm, a.{t) decreases 
to finite values and the density increases. However a(t) never reaches zero and 
attains a minimum value at t = tc given by 

tc = tm+ ^ ^ ^ . (62) 

For t > tc, a{t) starts to increase. We also note from eq.(l40l) that the curvature 
scalar i? — )■ — 00 as t — )■ and then starts increasing through finite negative 
values as t grows beyond tm- The singular behavior of the scale factor noted 
here should be contrasted with that of the FRW cosmological models in four- 
dimensional general relativity where the scale factor and energy density go to 
zero and infinity respectively as the initial moment is approached. 



5 Inflation 

The horizon problem in four-dimensional standard FRW cosmology is a conse- 
quence of deceleration in the expansion of the universe. The problem can be 
solved by postulating a phase of the universe, prior to the decelerating phase, in 
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which the expansion is accelerating and such a phase is called a period of infla- 
tion. Hence inflation is characterized by the following property for the scale factor 

a{t) : 

a(t) > . (63) 

Now as evident from the analysis of sect. 4, a > is readily achieved in our 
f{R) theory in two-dimensional spacetime and the universe is accelerating. The 
solutions obtained for the scale factor displayed power dependence on time akin 
to that of power-law inflation. It would seem that there is no need to require an in- 
flationary phase since matter or radiation dominated epochs yield an accelerating 
universe. Here we are not seeking to introduce scalar fields to propel acceleration 
as in the usual inflationary cosmology. We recall that one of the motivations for 
introducing modified or f{R) theories of gravity in four-dimensional spacetime is 
the desire to explain acceleration of the universe as an alternative to using scalar 
fields. For this purpose solutions for the cosmological field equations are sought 
in the absence of the matter fluid [5]. We carry out such an analysis in our case 
by considering solutions to eqs.(l32l) and (l33l) of sect. 3 with the R.H.S set equal 
to zero. We have earlier obtained a general solution for eq.(l33l) given by eq.dJTl) 
of sect. 3. However the parametric nature of that solution makes it difficult to use 
in eq.(l32l) in order to solve for a{t). Putting Ci = enables the integration in 
eq.(l37l) to be performed and leads to the solution given in eq. (|40l) which we write 
below as 

4n(2n - 1) 

^ - ~ (t-tp ' ^^^^ 

where n 7^ | , 1 and we have now denoted the integration constant by i. Using 
eq.(l64l) in eq.(l32l) with the R.H.S. set equal to zero yields the equation 

2n — 1 

a H ^ a = , (65) 

t-t 

the solution of which reads 

a{t) = A{t - ^ (66) 

where A > is a constant. We take the solution to hold for t > i . The cosmic 
acceleration is given by 

a(t) = 2n{2n - 1) A (t - t)"^"-! . (67) 
The Hubble parameter is 

a 1 — 2n 

H = - = , (68) 

a t — t 
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and 

For n < — 1 we can identify t with the onset of inflation t = U. Eq.(l66l) then 
describes a universe that expands with positive acceleration for t > ti . We also 
have H > and H < for t > which characterizes standard inflation. How- 
ever if we make the identification t = tj for n > 1, we will have a situation in 
which a(t) — t- oo as t — )• thus obtaining a universe that starts off already with 
an infinite size at the onset of inflation collapsing subsequently for t > U at an 
accelerated rate. Such a scenario can be avoided if i is instead taken to have a 
relatively large value so that t < i during the inflationary epoch. We write a{t) 
now as 

a(t) = A . (70) 

The universe then starts off with a relatively small non-zero size at t = U and 
expands with positive acceleration as time progresses. We also have 

2n-l 

H = ^ , (71) 

t-t 

1 - 2n 

so that H > and H < and we again have standard inflation. 

As in four spacetime dimensions we define the so called slow-roll parameter e by 

[5] 

e = -H/H'^ . (73) 

and in terms of which one has 

- = H^ + H=(l-e)H^ . (74) 
a 

Inflation can thus be attained only if £ < 1. In our present context e is given by 

1 



\2n - 11 ' 



(75) 



for both cases of n < —1 and n > 1 we clearly have e < 1. The slow -roll 
approximation corresponding to e « 1, then obtains when \2n — 1| >> 1. 

As we have stated above the solution for R given in eq.(l64l) arises as a special 
case of the general solution given in eq. (|37l) . As an alternative to solving eqs. (|32l) 
and (1331) one can derive an equation for the Hubble parameter [15], [5]. We write 
eq.(l32l)with the R.H.S set equal to zero 

2nRd + R'^a = 0. (76) 
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Now from eq.(l26l) of sect. 2 we obtain 



• 2 a 2d d ^^^^ 
R= + ^ ■ (77) 



a 



Substituting eqs.(l26l) and (1771) in eq.(r76l) one obtains 



n a d'd + a d^ + na^d = . (78) 



Next in terms of H, H and H we can express eq.(l78l). after some manipulations, 
as 

- nHH - 2{n - l)HH^ + + = . (79) 

It is customary, in dealing with equations such as this, to invoke the slow-roll 
approximation \H/H'^\ « 1 and \H/HH\ « 1, [15], [5]. Applying this to 
eq.dTOl) we obtain that 

-2{n~l)H + = . (80) 



The solution of eq.dSOl) is 



where t' is a constant. Eq.dSTl) in turn gives 

a(t) = A7(t-t')'("-i) . (82) 

with A' being another constant. Eq.(f82l) for a{t) is similar in structure to eq.dTOl) 
and the properties of the solution are therefore similar to what we discussed before 
and hence will not be considered any further. 

We shall next seek a general solution to eq.(l33l) for R{t) that holds for t close 
to the instant U that signifies the onset of inflation. Specifically we assume that 
t = ti is a regular point of eq. (l33l) and seek a solution for R{t) in the form of 
a power series confining ourselves to small values of t — ti. For simplicity we 
consider the case n = 2 for which eq. (I33l) becomes 

4i? + = . (83) 

We write 

oo 

R{t) = Y, bUt-Ur- (84) 

m=0 

Substituting eq. (|84|) in eq.(f83l) and solving we obtain 

b2 = -Ibl , (85) 
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^bo h , (86) 



etc. This leads to 



a{U) = Cexp{ - ^] ■ (90) 



R{t) = bo + h{t~ ti) - ^blit - Uf - ^bohit - U)' + (87) 

We remark that if inflation lasts for a short period of time then it is sensible to 
have a representation for R(t) as given in eq.(f87l). Moreover for sufficiently small 
t — ti we can approximate R(t) by the first two terms and substitute in eq.(|76l) 
with n = 2. Solving the resulting equation we obtain 

a(t) ~ C expj - ^[6o + h{t - U)f^ , (88) 
where C > is a constant. We can write eq.(l88l) as 

ait) ^ ai expj - ^ {[bo + bi{t - U)]^ - &^)| , (89) 

where 

n, = n.(t;\ = a PT/nx — 

I2bl^ 

From eq. ([88]) we obtain 

ix{t) = -^[bo + hi{t-U)fa{t), (91) 

and 

a(t) = I - \[ho + hi{t- ti)] + J^[6o + - U)f'^a{t) . (92) 

From eq. (f89l) we see that we must have 6i < to ensure that d > 0. We must also 
require a{t) to be increasing for t > t^. This can be achieved by having 6o > 
for then b^ + — tj) will start off at the value 6o and decreases reaching zero 
at t* — ti = — ^. During the interval, ti < t < t*, a(t) will be increasing. We 
must also require the cosmic acceleration d{t) to be positive during this interval 
and this leads to the condition 

^[6o + &i(t-t.)]'>l. (93) 

This inequality will continue to hold until t = tj < t* when a{tf) = 0. This 
implies that 

^[6o + 6i(t/-t,)]=' = l, (94) 
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which yields 

tf = U + bo\h\-'-2\b,\-'/\ (95) 

The time tj then signifies the end of inflation. Since R{ti) = bo and R{ti) = bi, 
the conditions 60 > and 61 < can be expressed as 

R{k) > (96) 

Riti) < . (97) 
We can also express the duration of inflation as 

tf-U = R{U)\R{ti)\-' - 2\R{U)\-'/' . (98) 

The Hubble parameter is given by 

H = -^h + bi{t-h)]' . (99) 

It thus decreases from an initial value Hi given by 

to a value Hf at the end of inflation where 

Hf = Hitf) = \b,\'/^ = \RiU)\'/''. (101) 

We note that 

H = -l[bo + b^it-U)], (102) 

is negative during U < t < tj and we thus have standard inflation. The slow-roll 
parameter is given by 

e = 8bl[bo + biit-U)]-'' . (103) 

We recall that for inflation to proceed one must have e < 1 and this leads to 
precisely the condition expressed in eq.(|93l) stated earlier. 

The number of e-foldings from t = ti to t = tj is defined by [5], [16] 



N= Hdt, (104) 



which is evaluated to give 

N 



(105) 
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In four dimensions, the solution of the horizon and flatness problems of big bang 
cosmology requires that iV > 70 , [17], [5]. If we assume that we can use this 
value in our two-dimensional universe, we find that 




(106) 



i.e., the Hubble parameter decreases to about 4.5% of its initial value by the time 
inflation ends. 

6 Quantization 

As we stated in the introduction two-dimensional spacetime models of gravity 
provide an arena where issues like quantization are studied since in such a set- 
ting they prove to be more tractable than in four-dimensional spacetime. In this 
section we thus consider quantization of the f{R) gravity theory defined by the 
action of eq.(3). Our objective is to derive the Wheeler-DeWitt equation for the 
wave function of the universe and obtain its solutions. Since we are considering 
a spatially homogeneous and isotropic universe, we drop the spatial intergral and 
write the action as 



We take for f(R) the expression given in eq.(31) and put n = 2. We use eq.(26) 
that expresses the scalar curvature in terms of the scale factor and write 



We notice the appearance of the second derivative of a in eq.(108). The stan- 
dard approach is to express the wave function in terms of a and it! [19]. Hence 
integrating by parts in eq.(108), we obtain 




(107) 




(108) 




(109) 



where 



^9m 

The canonical momenta are defined in the usual way 




(110) 




(111) 



da 



9. 
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R 



dR 



a 



9^ 



(112) 



The Hamiltonian is then obtained as 



Pah + PrR-^ 
1 



a 2g 



-aR . 



(113) 



Replacing Pa and Pr by and —i^ respectively in the Hamiltonian we obtain 
the Wheeler-DeWitt equation for the wave function of the universe 



9n 



— aR)ip(a,R) = . 



a dR da 2gf^ 
Instead of a and R we shall work with the variables 

^ = R + a, Tj = R — a . 

In terms of ^ and 77 the Wheeler-DeWitt equation becomes 



We seek solutions of eq.(l 16) in factorizable form 

and obtain the following equations for the functions X and Y 



(fx 



+ 



Q; ^2 q; 



^9l 

(PY 



ex^ 



+ 



9n 



-V' Y = Y , 



(114) 



(115) 



(116) 



(117) 



(118) 



(119) 



dv^ ^9l 9n 
where C is the separtation constant. The two equations are identical and hence it 
is enough to consider one of them. We first take a > and define 



7 



E 



a 



84 
Ca 
"2^ 



(120) 
(121) 
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In terms of 7^ and £■ eq.(l 18) reads 



^ + -I'^eX + 2 £; X = . (122) 

It is interesting to note that eq.(122) is identical to that describing the inverted or 
reversed oscillator discussed by several authors in a number of contexts [20] -[24]. 
By performing the change of variable 



y = V27e, (123) 

we cast eq.(122) into the form 

'^'^ + ^y^X + £X = , (124) 



dy'^ 4 

where e = Ej^. Eq.(124) is one of the standard forms of the equation for the 
parabolic cylinder functions. Two linearly independent solutions are given by the 
real functions W{e,y) and {e,—y) [25]. For \y\ » 1 and \y\ » \e\ these 
solutions display the following asymptotic behaviour 



2k /l 1 1 \ 

W{£, y 00) ^ J — cos y-y"^ + elny+-TT + -(pj , (125) 
1/ y 



W{e, y ^ -00) ~ y ^ sin (^^y^ + e ln\y\ + + ^(j)) , (126) 
where 

A;= (l + e-2"")^/2 -e""" , (127) 

<P^argr(^-ie) . (128) 

The functions W{e,y) and W{e, —y) satisfy the following normalization condi- 
tions [22] 



/ W{e,y)W{e',-y)dy^{ ^ Ji , (129) 

J -00 y (i+e-2T£)i/2 11 fc — fc . 

/oo 
W{e,y)W{e! ,y)dy = 2'K{\^e-'^'''fl'^b{e-e') . (130) 
-00 
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The parablolic cylinder functions can be expressed in several forms [25] and we 
can use the various relations between these forms to express W{a, x) in terms of 
the more familiar function Dp{x) for some p. In fact one can easily derive that 



where 



(131) 



^ = If-l7re + ^7r + i0) . (132) 



2 V 2 4 

Next we observe that the solutions to eq.(119) are identical to those of eq.(118) 
but expressed in terms of the vriable 77. Hence we can write the following for the 
wavefunction ^. 

V'(e,^) = V'i(OV'2(77) (133) 

where 

V'i(0=C'iiy(£;/7,v^0 + C2W^(^/7,- a/27 , (134) 



Mi) = C[ W{Eh, v/27 77) + C'^W{Eh, - v/27 77) . (135) 

We now consider the case in which the parameter a is negtive and write eqs. (118) 
and (119) as 

d'^X \a\ 2 C\a\ 

'"lyy^-^y. (137) 



We define 



^=-^, (139) 

and thus they retain the same forms as in eqs. (120) and (121) respectively. Fo- 
cussing on eq.(136) we write it as 

^^-^^eX^2EX . (140) 

In terms of y = ^, eq.(140) becomes 

f + (. + i-i.=)x^O, (141) 
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where 

1 E 

(7+- = . (142) 

2 7 

Eq.(141) has the form of Weber's equation [26] and possesses the following solu- 
tion 

X,{y) = DM = 2t+i y-^ (l/72) ■ (143) 

In the above W^,,^ is the Whittaker function. Expressing W^^i, in terms of the 
confluent hypergeometric function we can write 



r(i - £) 

For the second solution of eq.(141) we note that from eq.(143) giving the relation- 
ship betwen D^r and the Whittaker function, we know that D_(^_i(±iy) are solu- 

2 

tions linearly independent of D^{y) as W_a_i ^) is linearly independent 

2 2 4 4 

of Wo__^i_i (^). From the asymptotic behaviour of the confluent hypergeometric 
function we deduce that as y — > oo. 



X,{y)^e--y'' . (145) 
i.e., Xi 0. For y — > — oo we have 

My) - -^TT^e'-^V- y-''-' . (146) 
1 {-a) 

i.e., — > oo unless cr is a positive integer or zero in which case the R.H.S of 
eq.(146) vanishes. In fact we have the relationship 

D„(y) = 2-§e-^i/„(-^) , n=0,l,2, (147) 

that expresses the parabolic cylinder functions in terms of the Hermite poly- 
nomials Hn- Going back to eqs.(138), (139) and (142) with a — nwe obtain 

Cn = -^r-^ , (148) 

V2|«|2 

as the value of the separation constant. The functions that solve eq.(140) 

are then precisely those that describe the one-dimensional quantum oscillator. We 
write 
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The solutions Yn{ri) are identical in form and we obtain for the normalized wave 
function the following 

M^:V) = (^)' ^ ^"(^70 H^iV^v) e-i^«^+''^) . (150) 

For (T 7^ n the wavefunctions will not have finite norm and solutions of eqs.(136) 

and (137) of the type given in eq.(144) would have to be superposed, just as wave 
packets are constructed in quantum mechanics, in order to obtain wave functions 
capable of describing physical states. 
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7 Conclusions 



In this work we studied f{R) theories of gravity in two-dimensional spacetime 
with focus on applications to cosmology. With the metric taken to have to the 
FRW form we were able to obtain solutions for the cosmological field equations 
in the case of pure matter or radiation dominated universe when f{R) = R+aE!\ 
The remarkable feature of these solutions is that they readily describe an accel- 
erating universe in contrast to the standard FRW cosmology of four-dimensional 
general relativity. The horizon problem is also readily solved. As we have stated 
in sect. 2, the time evolution of the scale factor is not affected by the value of the 
curvature constant k. We have also seen that the solution for the radiation domi- 
nated universe, and one solution for the case of pure matter domination, describe 
a hot big bang. However an interesting solution in a matter dominated universe, 
given in eq.(l44l). describes a universe that kicks off with an infinite size and zero 
temperature at the start of matter dominance. It subsequently collapses to a finite 
size and then begins to expand. 

Now as we mentioned before, the interest behind the pursuit of f{R) theories is 
partially due to the desire to obtain a description of inflation without the intro- 
duction of scalar fields. This is done by seeking solutions to the cosmological 
field equations with the energy-momentum tensor set equal to zero, [5]. In sect. 
5 we obtained such solutions that characterize power law inflation. Furthermore, 
with inflation presumed to last for a short period of time, we obtained for the case 
n = 2 a solution for t near tj, the instant of onset of inflation. This solution dis- 
played exponential dependence on time. For this case we computed the duration 
of inflation and the number of e-foldings as well as an estimate for the change that 
ensues in the value of the Hubble parameter from the start to the end of inflation. 
The basic distinguishing feature between power law and exponential inflation ap- 
pears to be in the behavior of the Ricci scalar. Exponential inflation obtains when 
we assumed that R(t) can be expanded in a power series about t = U with finite 
coefficients. In particular R{ti) and R{ti) are finite. On the other hand in the case 
of power law inflation these quantities exhibit singular behavior att = i. Another 
characteristic of our inflationary solutions is that they do not depend on the pa- 
rameter '1 that appears in eq.(l3TI) for f{R). This is in contrast to the inflationary 
solution in four-dimensional f{R) theories where n = 2 describes the Starobin- 
sky model [18]. In that case with a written as a = g-g^, where the constant M 
has the dimension of mass, exponential inflation is obtained with a, H and R all 
depending on M, [5]. 

Interest in two-dimensional theories stems partially from the desire to investigate 
the quantum theory in a simple setting. Hence we carried out quantization of the 
theroy in the case of n = 2. The Wheeler-DeWitt equation was derived and its 
solutions obtained. We were able to solve the equation exactly in the entire do- 
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main of the variables, unlike the situation in the four-dimensional case [19], [27]. 
Interestingly we found that for a > the equation for the wave function coincided 
with that of the inverted oscillator. For a < the wave function, under certain 
conditions, turned out to be a product of two quantum harmonic oscillator wave 
functions in the variables ^ = i? + a and r] = R — a. In conclusion we have stud- 
ied some aspects of classical and quantum cosmology in two-dimensional f{R) 
theories. Clearly a lot more issues need to be investigated and we hope to return 
to them in the near future. 
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